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In this paper, we briefly review the concept of memory circuit elements, namely memristors, 
memcapacitors and meminductors, and then discuss some applications by focusing mainly on the 
first class. We present several examples, their modeling and applications ranging from analog 
programming to biological systems. Since the phenomena associated with memory are ubiquitous 
at the nanoscale, we expect the interest in these circuit elements to increase in coming years. 



I. INTRODUCTION 

Currently, there is strong interest in circuit elements 
with history-dependent properties. These form a general 
class of memory elements^ which includes memory resis- 
tors (memristors)^, memory capacitors (memcapacitors) 
and memory inductors (meminductors). Quite gener- 
ally, memristors, memcapacitors and meminductors can 
be considered as time- and history-dependent generaliza- 
tions of their standard counterparts. 

Mathematically, an nth-order li-controlled memory el- 
ement is defined by the equations^ 



y(t) = g{x,u,t)u{t) 
X = f {x^ t) . 



(1) 

(2) 



Here, u{t) and y{t) are any two circuit variables (current, 
charge, voltage, or flux) denoting input and output of 
the system, x is an n-dimensional vector of internal state 
variables, ^ is a generalized response, and / is a contin- 
uous n-dimensional vector function. Special interest is 
devoted to devices determined by three pairs of circuit 
variables: current- voltage (memristors), charge- voltage 
(memcapacitors), and flux- current (meminductors). Two 
other pairs (charge-current and voltage-flux) are linked 
through equations of electrodynamics, and therefore do 
not give rise to any new device. Devices defined by the 
relation of charge and flux (which is the integral of the 
voltage) are not considered as a separate group since such 
devices can be redefined in the current- voltage basis^. 

We stress first that these memory devices are dynam- 
ical systems, namely systems whose state and response 
may change in time. In typical applications, memory 
devices operate under the action of a time-dependent 
(not static) input. It is also worth noting that the re- 
lations ([1]) and ([2|) are quite general. In particular, they 
may define systems whose properties cannot be repro- 
duced with standard resistors, capacitors and inductors. 
In other words, there is no possible combination of stan- 
dard circuit elements that can reproduce the dynamical 
properties of certain memory elements. In this sense, 
memristors, memcapacitors and meminductors can be 
considered as "fundamental" circuit elements. However, 
at least in the case of memcapacitors and meminduc- 
tors, the relations ^ and ([2]) may also represent some 



other systems whose dynamical properties can be simu- 
lated by an appropriate combination of standard (albeit 
possibly non-lineax) elements (for examples of these, see, 
e.g., Refs. i andy). Nonetheless, these latter memory 
systems are still of great importance since they provide 
a complex functionality within a single electronic struc- 
ture^. 

The area of research of memristors is more advanced 
both technologically and theoretically than that of mem- 
capacitors and meminductors. In fact, although some 
experimental systems have been identified as memcapac- 
itiveii^r^ and meminductive^, their number is still small 
and possible applications for them are less developed. In 
this review, we will then focus on memristors^. 

It is convenient to introduce voltage-controlled and 
current-controlled memristors^. An nth-order current- 
controlled memristive system is described by the equa- 
tions 



VM{t) = R{x,I,t)I{t) 
X = f{x,I,t) 



(3) 

(4) 



with X again a vector representing n internal state vari- 
ables, Vm(^) and I{t) denote the voltage and current 
across the device, and is a scalar, called the memris- 
tance (for memory resistance). The equation for a charge- 
controlled (or ideal) memristor is a particular case of Eqs. 
([3]) and when R depends only on charge, namely 



VM=R{q{t))I, 



(5) 



with the charge related to the current via time derivative: 
/ = dq/dt. 

We can also define an nth-order voltag e-contiolled 
memristive system from the relations 



I{t)=G{x,VM,t)VM{t) 
X = f{x,VM,t) 



(6) 
(7) 



where we call G the memductance (for memory conduc- 
tance). 

There are several common properties of memris- 
torsii^i^. They are passive devices (provided R > or 
G > 0) without an energy storage capability (unlike 
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FIG. 1. Schematic of typical I — V curves of a memristor 
obtained when an ac- volt age of frequency / is applied. At 
low frequencies fiow^ the I — V curve is a non- linear curve, 
at intermediate frequencies flow < f < fhigh, it is a pinched 
hysteresis loop, and at high frequencies fhigh, it is a straight 
line. Note that all curves pass through (0,0). 



memcapacitors and meminductors which can store en- 
ergy). Typically, at high frequencies, they behave as lin- 
ear resistors and at low frequencies as non-linear resistors 
(assuming existence of a steady-state solution of Eq. dH 
or ([7j)). A distinctive feature of well-defined memristors 
(with non-zero R or G) is a pinched hysteresis loop on 
the current- voltage plane (see Fig. [1]). 

Because of the short format of this review, we are not 
able to include all topics concerning memristors or to 
cover the ones included in due detail. We will focus 
mainly on our results in this field and organize the text 
in the following order. Sec. |TT1 describes several differ- 
ent realizations of memristors. A theoretical modeling 
is discussed in Sec. Illli We overview several important 
applications of memristors in Sec. IIVI Finally, Sec. |Vl 
presents our conclusions. 



II. EXPERIMENTAL REALIZATIONS 

Many different systems exhibit memristive behavior. 
Early-identified memristor's examples^ include thermis- 
tors^ and ionic systemsiS. Recent interest in mem- 
ristor's realizations has centered on nanoscale systems 
and nanostructures as we discuss below. Some of them 
were intensively studied in the context of development of 
resistive- switching memory during last 10-15 years. This 
area, to the best of our knowledge, was initiated by Hick- 
mott in 1962 by the observation of hysteretic behavior in 
oxide insulators^^. 

At the present time, it is well known that memristive 
effects can be observed in such diverse classes of materials 
as binary oxides (Ti02, CuO, NiO, CoO, FesOa, MoO, 
V02)^'^^^^, perovskite-type oxides (Pri-a^Caa^MnOs, 
SrTiOsiCr)!^ sulfides (Cu2S,Ag2S)Sa-24^ semiconduc- 
tors (Si, GaAs, ZnSe-Ge)^^^'^, spintronics materials^^^^^ 



and organics^— . In this review, we limit our discus- 
sion to the binary oxide Ti02 and spintronics systems. 
We will also discuss a memristor emulator—*^ - an elec- 
tronic module simulating the memristor behavior. 



A. Binary oxides materials 

Ti02 is an example of binary oxide showing bipolar 
resistive switching. Recently, a memristive mechanism 
of its switching behavior has been identified^^^^^. It was 
suggested that the switching involves changes to the elec- 
tronic barrier at the Pt/Ti02 interface induced by drift 
of oxygen vacancies under an applied electric field. When 
vacancies drift toward the interface, they create con- 
ducting channels that short-circuit the electronic barrier. 
When vacancies drift away from the interface they elimi- 
nate these channels, restoring the original electronic bar- 
rier This material demonstrates promising character- 
istics for future ultra-high density memory applications 
such as fast read/write times (~10ns), high on/off ratios 
(~ 10^), suitable range of programming voltages, and 
possibility to fabricate small size cells. 



B. Spintronic materials 

Semiconductor spintronic systems^^ were identified as 
memristive systems by Pershin and Di Ventra^^. In semi- 
conductor spintronics, memory effects are completely due 
to the electron spin degree of freedom. When an ex- 
ternal control parameter changes, it takes time for the 
electron spin polarization to adjust to a new control pa- 
rameter value (typically, the equilibration is governed by 
electron-spin diffusion and relaxation processes). There- 
fore, if a system's resistivity depends on the level of elec- 
tron spin polarization (as in semiconductor-ferromagnet 
junctions^^'^^) then such a system is fundamentally mem- 
ristive^^. Moreover, memory effects were predicted for 
spin Hall effect systems with inhomogeneous charge den- 
sity distribution transverse to the driving electric field^^. 
It is interesting that in this case memory effects can be 
observed directly in the transverse voltage. Fig. [2] shows 
transverse voltage hysteresis loops that demonstrate typ- 
ical memristive behavior: non-linear dependence at low 
frequencies, pronounced hysteresis at higher frequencies 
and hysteresis collapse at very high frequencies^^. 

In metal spintronics, spin torque transfer known for 
a long time^^^^ can also be categorized as memristive. 
In spin torque transfer systems, the resistance is deter- 
mined by the relative magnetization between opposite 
sides of a magnetic tunnel junction. Current flowing 
through the junction induces spin torque, in turn chang- 
ing the relative magnetization. Krzysteczko et al.^^ fab- 
ricated a modified magnetic tunnel junction with a re- 
sistive switching material inside. A spintronic memristor 
based on spin-torque-induced magnetization motion was 
discussed by Wang et alM. 
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demonstrated^^. 



III. THEORETICAL DESCRIPTION 

For practical use of memory elements, the ability to 
predict their dynamics in electronic circuits is important. 
For this purpose, it is necessary to know the functions 
gix^u^t) and f{x^u^t) entering Eqs. (jl|2p that can be 
later used, for example, in SPICE modeling^^^^. Gen- 
erally, these functions can be constructed phenomeno- 
logically or derived from an appropriate physical analy- 
sis. Although models of resistive- switching devices were 
suggested before 2008"^-, the realization that these de- 
vices belong to the class of memristive systems^^ has 
spurred further progress in understanding their behav- 
ior—"—. Below, we consider a phenomenological model 
of memristor suggested by two of us (YVP and MD) to 
explain the learning properties of unicellular organisms^^- 
and a computationally simple tunneling barrier model. 
The tunneling barrier model includes both activation- 
type of internal state variables dynamics and electron 
tunneling through the barrier and it is indeed a physical 
realization of the above phenomenological model. It is 
also useful to understand the memory resistive proper- 
ties of several systems of present interest. 



FIG. 2. Spin memristive effects in a semiconducting system 
with inhomogeneous electron density in the direction perpen- 
dicular to main current flow^^. Here, we show the transverse 
voltage as a function of applied electric field at different ap- 
plied field frequencies as indicated. (From Ref. [2^ with per- 
mission. Copyright American Physical Society, 2009). 



C. Memristor emulator 

A memristor emulator is an electronic module com- 
posed of standard (off-the-shelf) electronic components 
whose circuit response imitates that of memristor. Mem- 
ristor emulator consisting of a microcontroller-controlled 
digital potentiometer and an analog-to-digital converter 
was recently demonstratedS^^. Using the analog-to- 
digital converter, the microcontroller obtains informa- 
tion about the voltage applied to the memristor emulator 
and updates the digital potentiometer resistance using 
equations of voltage-controlled (Eqs. ()6|7p ) or current- 
controlled (Eqs. (|3|4p ) memristive systems. From the 
point of view of the external circuit, the memristor em- 
ulator is then a black box with two connections that 
operates exactly as a memristor. Its low cost, simplic- 
ity and flexibility make it useful for building prototype 
circuits with memrsistors for both research and educa- 
tion. Moreover, using the memristor emulator (which 
could be replaced by a real memristor), modules showing 
memcapacitive or meminductive behavior have also been 



A. Phenomenological model 

A phenomenological model of a memristor has been 
suggested in Ref. [53|. This model captures two very im- 
portant experimental features of Ti02 memristors^^: i) 
the fact that the resistance changes between two limiting 
values Rmin and Rmaxi and ii) a threshold- type resis- 
tance dynamics. The latter is manifested in the following 
experimental observatioi>i^: at low applied voltages the 
resistance of the memristor changes slowly while at high 
applied voltages its change is fast. The phenomenolog- 
ical model describes a voltage-controlled memristor (see 
Eqs. dSlIll)) with Rm = x m Eq. ^ and Eq. ^ written 
as (with the memristance acquiring the limiting values 

Rmin and Rmaxy^ 



X = WVm + 0.5 {a - p) [\Vm + Vt\ - \Vm - Vt\]) 

X0{X- Rmin) {Rmax - x) , (8) 

where a and P are constants defining the memristance 
rate of change below and above the threshold voltage Vr; 
Vm is the voltage across the memristor, and 6{-) is the 
step function. This model has been used in several cases 
ranging from modeling of learning to emulators^^^^^^^^. 
Moreover, similar models can be applied to memcapaci- 
tors and meminductorsi. 
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FIG. 3. Schematic of Ti02 memristor (a) and of its resistivity switching mechanism (b). The Schottky barrier is reduced in 
reverse bias situation when oxygen vacancies are shifted to the left. 



B. Tunneling Barrier Model 

Quite generally, let us assume that the electron trans- 
port through a voltage-controlled memristor involves a 
quantum mechanical tunneling through a barrier whose 
width and height are determined by internal state vari- 
ables of the device. A particular example of application 
of this model regards the Pt/Ti02/Pt thin- film memris- 
tor^^, in which a Schottky barrier plays the role of the 
tunneling barrier and the internal state variable is the 
density of oxygen vacancies in the Schottky barrier re- 
gion. For the sake of clarity we then focus on this par- 
ticular system. We also note that the Schottky barrier 
was not exphcitly included in previous titanium dioxide 
memristor models^'—. 



1. Dynamics of Internal State Variables 

The structure of a Ti02 memristor is very simple^^, 
as shown in Fig. [3fa). It consists of two metallic (e.g., 
platinum) electrodes separated by a titanium dioxide thin 
film containing a region with a high concentration of oxy- 
gen vacancies, that we denote with Ti02-x, with x de- 
scribing the vacancies concentration. The right contact 
in Fig. [3l^a) of the Ti02-a; region with Pt can be consid- 
ered as ohmic^^. The left contact (Ti02 with Pt) involves 
a Schottky barrier 

A fraction of oxygen vacancies can be moved from the 
TiOo -X region to the rest of Ti02 by an applied elec- 
tric field via the phenomenon of electro-migration (see, 
e.g., Ref. |53). Let us consider only mobile vacancies 
and track only the total number of mobile vacancies in 
the Schottky barrier region that extends inside of Ti02. 
Moreover, we assume that in our device the separation 
between Pt electrodes is of the order of two Schottky bar- 
rier lengths (here, by Schottky barrier length we mean 
its length when the device is in its high-resistance state). 
Only the vacancies in the Schottky barrier region are im- 
portant to us since they modify the shape of the barrier 



(this is schematically shown in Fig. [3jb)) influencing di- 
rectly the electron transport through the device^-. 

The vacancy dynamics can be described by a rate equa- 
tion 



dNL 
dt 



1 

r 



-e{y)NL + e{-v)NR] , 



(9) 



where V is the applied voltage and Nl^r are the number 
of mobile vacancies in the left L (Schottky barrier) region 
and right, i?, remaining part of titanium dioxide. The 
total number N of mobile vacancies is assumed to be 
constant {N = Nl-\-Nji). At the initial moment of time, 
Nl can take any value between and N depending on 
the initial state of the system. Limiting cases Nl = and 
Nl = N correspond to high-resistance and low-resistance 
states of the memristor, respectively. 

The time r taken for a vacancy to move from one re- 
gion to another can be found using an activation-type 
expression for the drift velocity as proposed in Ref. [50| 
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\v\ 2/j.Eo sinh(|£;|/£;o) 2/j.Eo sinh (| £;o) 

/i Eo 



exp 



\LBa) 



if \V\ > LEo 



(10) 



where /i is the ion mobility, = ^ksT / {qa)^ T is the 
temperature, q is the elementary charge, a is the poten- 
tial period, and for the electric field E we have used the 
approximation |£^| = |V^|/i^, where L is the length of ti- 
tanium dioxide material. The derivation of Eq. (p!Q|) is 
based on the assumption that the oxygen vacancies are 
confined by a periodic potential due to all other atoms 
in the solid along the migration path. The application of 
the external electric field deforms the potential reducing 
the potential barrier between adjacent potential wells^^. 
As a result, the expression for the drift velocity becomes 
very non- linear at strong applied fields^^. Using Eq. (fTOj) , 
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Eq. (j9j) can be rewritten as 



dt 



L 



smh[\V\/{LEo)] 



-0{V)nL ^ 0{-V) {1 - ul}] , 
(11) 



where ul = Nl/N. A numerical solution of Eq. (fTTj) 
is shown in Fig. 3^. Later, we use nL{t) as an input 
parameter for tunneling current calculations. 



2. Electron transport 

Experimental data^^ suggest that the tunneling cur- 
rent through the Schottky barrier is the main limiting 
factor to electron transport. This current can be approx- 
imated by a standard expression^ 



Ae-^^ ircikT 
{cikTf sin(7rci/cT) 



[1 



-ciqVl 



(12) 



where A = ATrm^ q{kBT)'^ / is the classical Richardson 
constant, m* is the effective electron mass, bi and ci are 
constants. The expressions for bi and Ci are different for 
forward and reverse biases. In particular, in the forward 
bias when electrons tunnel from Ti02 into the metal (ex- 
pressions for bi and ci for the case of reverse bias can be 
found in Ref. 56), 



bi 



EB-qV 



Cl 



2-Eoo 



In 



A{Eb - qV) 



(13) 
(14) 



where Eb is the barrier height, difference between work 
function of the metal and the affinity of the semicon- 
ductor, ^ is the Fermi energy of the semiconductor and 
the energy ^oo is a material-dependent parameter that 
is inversely proportional to the Schottky barrier width^^. 
For the sake of simplicity, a constant barrier height is 
assumed in our calculations. The tunneling through the 
barrier dominates thermionic emission and thermionic- 
field emission mechanisms when Eqq ^ kT (the regime 
we consider here)^^. 

We assume that the Schottky barrier width, changes 
linearly with the number of the oxygen vacancies drifted 
to the Schottky barrier region of the device, that is 
/ = /o(l — ariL) where < a < 1 is a constant. There- 
fore, the maximum width of the Schottky barrier is Iq 
and the minimum width is /o(l — <^)- As the energy Eqq 
is inversely proportional to / (see Ref. 56), we can write 
that £^00 = ^oo/(l — (^ul) where £^qq is the value of Eqq 
at = 0. Fig. [JJd depicts memristor's I — V curves ob- 
tained using Eq. ([12]). These curves have a shape typical 
of the memristor pinched hysteresis loop as anticipated 
in Fig. [Hand typical memristor's frequency behavior. 
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FIG. 4. Simulations of an AC-biased memristive device: (a) 
time dependence of the applied voltage V, vacancies concen- 
tration ul and current density j, and (b) I — V curves ob- 
tained when different frequency ac signals are applied to the 
device. These plots were obtained using the following set of 
parameters: V{t) = Vbsin(27r/t) with Vo = IV, / = lO'^Hz in 
(a) and as specified on the plot in (b), Eb = 1.5eV, T = 300K, 
^ = 0.2eV, E^o = 2eV, a = 0.8, L = lOnm, = 2 • lO^V/cm, 
^ = 3-10"^VV(Vs). 



IV. APPLICATIONS 



Potential applications of memristors are envisioned in 
both digital and analog domains. Of great importance 
among digital applications are non-volatile solid-state 
memory, signal processing and programmable logio^*^. 
Analog applications of memristors are based on the pos- 
sibility to continuously vary their resistance and, there- 
fore, on their ability to store more information than in 
the digital regime. Analog applications include analog 
signal processing, learning circuits, programmable ana- 
log circuits and neuromorphic circuits. Below we discuss 
several analog applications of memristors^ '^^^ . 
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FIG. 5. (a) Learning circuit^^ is composed of four circuit 
elements: resistor R, capacitor C, inductor L and memristor 
M. The input voltage V{t) corresponds to varying environ- 
mental conditions, the voltage on capacitor corresponds to 
the response signal. For more details, see Ref. 53. (Adapted 
from Ref. 53 with permission. Copyright American Physical 
Society, 2009) . (b) A hypothetical memristive neural network 
involving neurons (circles) connected by memristive synapses. 



A. Learning circuits 



Quite generally, a "learning circuit" is an electronic 
circuit whose response at a given time adapts according 
to signals applied to the circuit at previous moments of 
time^^. Memristors are ideal components for such a cir- 
cuit as they provide non- volatile information storage and 
compatibility (as a time-dependent resistor) with other 
circuit elements. Two of us (YVP and MD) have recently 
suggested^^ a learning circuit which mimics adaptive be- 
havior of the shme mold Physarum polycephalum from 
the group of amoebozoa. This work was inspired by an 
interesting experimental observation^^: when Physarum 
polycephalum is exposed to a pattern of periodic envi- 
ronmental changes, it learns and adapts its behavior in 
anticipation of the next stimulus to come. It was demon- 
strated^^ that such behavior can be described by the re- 
sponse of a simple electronic circuit shown in Fig. [5^ 
composed of an LC contour and a memristor in parallel 
with the capacitor. The memristive function employed 
in that work is Eq. [51 When a periodic signal is applied 
to the learning circuit, the voltage across the capacitor 
significantly changes and can exceed the threshold volt- 
age of the memristor. This leads to an increase in the 
resistance of memristor and, consequently, in a smaller 
damping of the LC contour. Therefore, the LC contour 
oscillations are maintained for a longer period of time in 
analogy with the same type of behavior as the amoeba's 
when subject to periodical environment changes. 



B. Neuromorphic circuits 

Neuromorphic circuits are circuits whose operation is 
meant to mimic that of the brain. In these circuits, mem- 
ristors can be used as synapses whose role is to provide 
connections between neurons and store information. The 
small size of solid-state memristors is highly beneficial for 
this application since memristors' density can be of the 
same order of magnitude as the density of synapses in 
human brains^^. Therefore, using memristors, the fabri- 
cation of an artificial neural network of a similar size of 
a biological brain becomes possible. 

Typically, an artificial neural network consists of many 
artificial neurons connected by artificial synapses as we 
show schematically in Fig. [SJd. The first memristor-based 
neural network realized experimentally^^ contained only 
three neurons connected by two synapses. However, we 
(YVP and MD) have shown that even such a small net- 
work could perform a fundamental property of the human 
brain: associative memory. Associative memory is closely 
related to the so-called Pavlovian training^^ in which a 
particular response to a given stimulus develops. The 
most notable experiment of associative memory is that 
of a dog to which food is shown and, at the same time, the 
tone of a bell is rang so that, with time, the dog salivates 
at the ring of the bell only. In the electronic 3-neuron 
network demonstrated in Ref.^^, two input neurons were 
responsible for the "sight of food" and "sound" events, 
while the output neuron generated a "salivation" com- 
mand. It was shown that, if one starts with an untrained 
state of the synapse connecting the input "sound" neu- 
ron and output neuron and exposes the neural network to 
"sight of food" and "sound" signals simultaneously then 
an association between "food" and "sound" develops and 
an output signal is generated when only the "sound" in- 
put signal is applied. 



C. Programmable analog circuits 

In programmable analog circuits, memristors can be 
used as digital potentiometers^. The main idea is to 
apply small amplitude voltages to memristors when they 
are used as analog circuit elements and high amplitude 
voltage pulses for the purpose of memristor's resistance 
programming. Since the state of memristor appreciably 
changes only when the voltage applied to it exceeds a cer- 
tain threshold^^, the resistance of memristor is constant 
in the analog mode of operation and changes by discrete 
values with each voltage pulse. Using this idea, several 
programmable analog circuits demonstrating memristor- 
based programming of threshold, gain and frequency 
were demonstrated 
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V. CONCLUSIONS AND OUTLOOK 



As this short review has shown, memory circuit el- 
ements are new promising components for future elec- 
tronics. Their advantage is based on a combination of 
a history-dependent behavior with properties of a basic 
circuit element (such as resistance, capacitance or induc- 
tance). Of equal importance is the possibility of realizing 
memory elements at the nanoscale. The cost of certain 
implementations of these elements should be low because 



of their simple structure. We therefore anticipate that as 
memristors, memcapacitors and meminductors possibly 
develop into commercially available products in the com- 
ing years, novel exciting applications will be developed 
with impact not only for technology but also for funda- 
mental science. 
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